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Abstract 

We investigate numerically on the lattice the interplay of universality classes of the three- 
dimensional Yukawa model with U(l) chiral symmetry, using the Binder method of finite 
size scaling. At zero Yukawa coupling the scaling related to the magnetic Wilson-Fisher 
fixed point is confirmed. At sufficiently strong Yukawa coupling the dominance of the chiral 
fixed point associated with the 3D Gross-Neveu model is observed for various values of the 
coupling parameters, including infinite scalar selfcoupling. In both cases the Binder method 
works consistently in a broad range of lattice sizes. However, when the Yukawa coupling is 
decreased the finite size behavior gets complicated and the Binder method gives inconsistent 
results for different lattice sizes. This signals a cross-over between the universality classes of 
the two fixed points. 
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1 Introduction 



Some strongly coupled lattice field theories in 4 dimensions (4D) possess perturbatively 
unaccessible critical points where scaling properties are understood only poorly or not at 
all. Examples are noncompact QED [jl|, compact QED without matter fields (pure QED) 

or with fermions , gauged Nambu-Jona-Lasinio or Yukawa models and models with 
fermions, gauge field and charged scalar at strong gauge coupling A clarification of their 
critical behavior and of the continuum limit taken at such points is desirable at least for 
two reasons: Firstly, the fundamental question of the existence of 4D quantum field theories 
defined on nongaussian fixed points has never been settled. Secondly, finding a 4D theory 
interacting strongly at short distances could contribute to the development of theoretical 
scenarios for dynamical symmetry breaking as possible alternatives to the Higgs mechanism 
in the standard model and its extensions. 

Except the pure QED, a chiral phase transition, with the chiral condensate (xx) 
an order parameter, takes place at all the critical points mentioned above. But it always 
differs in some qualitative way from the classical model for chiral symmetry breaking, the 
Nambu-Jona-Lasinio model. This is encouraging, as that model is even nonperturbatively 
nonrenormalizable ^ and thus of very limited use. The differences consist mainly in an ad- 
mixture of some other phenomena like confinement, monopoles, magnetic or Higgs transition, 
additional states of vanishing mass, etc., intertwining with the chiral transition, but occur- 
ing also in other situations, including those without fermions. This increases the hope for a 
fundamental difference from the Nambu-Jona-Lasinio model, but also makes the transitions 
perplexingly complex and difficult to analyze. In particular, the genuine character of the 
transition might be hidden behind some prescaling phenomena caused by some component 
of the mixture, or by a crossover between different universality classes. 

In this paper we study the interplay of the chiral and magnetic phase transitions in a 3D 
lattice Yukawa model (Y3 model) with global U(l) chiral symmetry as an exercise for the 
investigations of analogous but more complex situations in 4D. The Y3 model has nontrivial 
fixed points, a property searched for in 4D. We would like to learn how to detect such points, 
and what are the possible obstacles when the scaling properties are investigated numerically 
in the situation of intertwining phenomena. 

The couplings of the Y3 model are the scalar hopping parameter k, the scalar quartic 
selfcoupling A and the Yukawa coupling y. The action is given in subsec. 2.1. The phase 
diagram is shown schematically in fig. |I|. We concentrate on the transition between the 
paramagnetic (PM) and ferromagnetic (FM) phases. The two-dimensional PM-FM sheet of 
2nd Qr(-[gr phase transitions connects the critical line of the purely scalar </)| model at y = 
and the critical point of the 3D Gross-Neveu (GN3) model at n = X = 0. On this sheet the 
Y3 model is expected to have two nontrivial fixed points: 

1. Wilson-Fisher fixed point (WFfp) [0 of the pure scalar 2-component 03 theory, whose 
most familiar representative is the 3D XY (XY3) model. The phase transition is of 
magnetic type. 

2. Chiral fixed point {xfy), most naturally associated with the GN3 model with U(l) 
global chiral symmetry and a chiral phase transition. The existence of this fixed point 
is related to the nonperturbative renormalizability of the GN3 model (see ||^ and 
references therein). 

The sketch of the renormalization group fiow in fig. |^ represents a plausible scenario for 
what happens along the critical PM-FM sheet: The magnetic WFfp describes only the 0| 
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Figure 1: Schematic phase diagram of the Y3 model. The region below the upper critical 
surface is the paramagnetic phase (PM), the region above it the ferromagnetic phase (FM). 
The y = plane and the /t = A = line correspond to the (/)| and GN3 models, respectively. 
For negative values of the parameter k we further expect an antiferromagnetic phase (AFM) 
and a ferrimagnetic phase (FI). We have investigated the PM-FM transition for ~ and 
K > 0, in particular at the points and directions indicated by the arrows. 



theory. The xfp presumably dominates (has a domain of attractivity) everywhere as long as 
the Yukawa coupling does not vanish, and in the limit of infinite cutoff the Y3 model is thus 
equivalent to the GN3 model. This expectation has been recently supported at weak scalar 
self coupling and large Yukawa coupling by the 1/N expansion ^ and a consequent 

combined analytic and numerical investigation A discussion of the equivalence between 
the Yukawa and four-fermion theories, as well as earlier references, can be found in ref . [|13| . 

In fig. ^ we show schematic RG flows also outside the critical sheet for three special cases 
of restricted parameter space: y = 0, k = and k = X = 0. This figure indicates that the 
known RG flows in the (p^ and GN3 models can be consistently embedded into the RG flows 
in the Y3 model. 

When in the Y3 model the Yukawa coupling decreases and the (j)^ theory is approached, 
the WFfp gets presumably influential, as some crossover to the magnetic universality class 
must occur. This consideration warns us that for limited lattice volumina and consequently 
limited correlation lengths either no unique finite size scaling behavior can be found or the 
wrong fixed point dominates. Thus a detection of the genuine - presumably chiral - character 
of the transition gets more and more difficult in numerical simulations. This is the situation 
we are most interested in, as it might occur in 4D without a prior warning. 

Apart from this particle physics motivation our work might be of interest also for other 
reasons, related mainly to statistical mechanics: 



1. We have applied the Binder method of finite size scaling analysis |15[ to the chiral 
phase transition and found that it works very well also when a composite scalar field 
is used in the finite size scaling theory, as long as the xfp alone dominates the finite 
size scaling behavior. 
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Figure 2: A suggestion for the renormalization group flow on the PM-FM critical surface of 
the Y3 model. The fixed points are Gfp (Gaussian), WFfp (Wilson-Fisher) and xfp (chiral, 
or GN3). The indicated position of the x^P is very schematic, it could lie anywhere on the 
PM-FM sheet, at y > 0. 




Figure 3: The schematic RG-flow in the (f)'^ model (y = 0, k > 0), in the k = surface of 
the Y3 model and on the k = A = line, which corresponds to the GN3 model. The fixed 
point of the latter model is indicated by an arrow. 



2. A transition between various universality classes in finite volumina has been investi- 
gated recently [0, in some spin models, but, to our knowledge, until now in no 



models with fermions. Thus we make a new contribution to the experience with this 
sort of complex finite size behavior. As in spin models, it is the failure of the Binder 
method which indicates a change of the universality class. 



3. Sometimes an intermediate universality class could exist |T^. This would be very 
surprising for the Y3 model, nevertheless we have verified that this is most probably 
not the case here. 

We now briefiy describe the contents of the paper and the main results: 
In the next section we introduce the Y3 model and determine its phase diagram (fig. |l|), 
both by means of the effective potential in the one loop approximation, and by performing 
numerical simulations on a small lattice at many points in the three-dimensional parameter 
space. The most useful order parameter is the scalar field expectation value, even if this 
field can be considered as composed of a fermion pair. We mention some results on the 
fermion and boson masses both in the symmetric phase and in the phase with broken chiral 
symmetry. 

In sec. 3 we shortly review the Binder method allowing a determination of several critical 
exponents by an analysis of finite size effects. The most useful exponent is the correlation 
length exponent v obtained from the Binder-Challa-Landau (BCL) [|I4| , |18] cumulant. 



The magnetic transition of the (\)\ theory is investigated in sec. 4. After localizing the 
critical line we concentrate on the case A = 00 (the XY3 model) and a case of an intermediate 
scalar selfcoupling (A = 0.5). The obtained exponents are consistent with each other and with 
the value expected from analytic investigations of the WFfp [v = 0.67). Also the values of 
the renormalized coupling extrapolated to infinite cutoff are consistent. The Binder method 
is compared with two other approaches to finite size scaling and found to be most suitable 
for our purposes. 

Sec. 5 deals with the chiral transition in the GN3 model at A = both in the auxiliary 
scalar field formulation (k = 0), and with a dynamical scalar field (k varied and y kept at the 
critical value, y = 1.09). In both approaches to the critical point the Binder method works 
comparably well for all the lattice sizes we used (6^ - 24^) and gives consistent results for 
critical exponents. In particular, u = 1.03(11), which is a value consistent with theoretical 
expectations |jl3|, |1^, |2^ and significantly different from the value found for the 03 model 



at ?/ = 0. Thus the difference between the magnetic and chiral universality classes is clearly 
observed in the y = and A = limit cases. Their common property is that the Binder 
method works in an exemplary way in the whole range of lattice sizes we used. 

In sec. 6, the Y3 model with a large Yukawa coupling, y = 1.1, is investigated at the 
maximal value of the scalar selfcoupling A = 00. Also here the Binder method works quite 
well, and we find u = 0.88(6), a value slightly lower than, but within errors still consistent 
with that found in the GN3 model. This confirms the appurtenance of the Y3 model with 
both couplings y and A strong to the same chiral universality class as the GN3 model, and 
thus the physical equivalence of both theories. 

However, difficulties arise when the Yukawa coupling decreases. As we describe in sec. 7, 
at A = 00 and y = 0.6 the BCL cumulants cross at different points when only small (6^ - 
10'^) or large (10'^ - 24^) lattices are considered, suggesting different values of the critical 
K. Restricting ourselves to the larger lattices only, we find the Binder method to work, 
giving i> = 0.99(23). This value is consistent with the GN3 model value, but has a large 
error. On small lattices the obtained value of u is significantly lower and close to the value 
in the 0| model. As we describe in detail in the same section, at A = 00 and y = 0.3 the 



Binder method gives inconsistent results in the whole range of lattice sizes 6'^ - 32'^ we have 
investigated. This can be interpreted as a situation in which none of the two fixed points 
alone dominates the finite size effects on lattices of these sizes, i.e. as an interplay of or a 
crossover between universality classes. We find no sign for the existence of an intermediate 
universality class. 

As we conclude in sec. 8, an interplay of magnetic and chiral phenomena in the Y3 model 
thus results in uncontrollable finite size effects. However, inconsistencies in the application of 
finite size methods become apparent only when a broader range of lattice sizes is investigated. 
This might serve as a warning for investigations of critical points with a mixture of chiral 
and some other critical behaviour in 4D lattice field theories. 



2 The Y3 model and its phase diagram 
2.1 The action 

In order to investigate the breakdown of a continuous chiral symmetry we use staggered 
fermions In the lattice parametrization the action of the Y3 model is 



S — 5*3 + Sp + 
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= §J2^iJ2i'f''^+b + ^ex<PU,)xi, (2.1) 



where the integer 3-vectors x, x + fi and x + b denote, respectively, lattice site, its nearest 
neighbors and corners of the associated elementary cube (both in positive direction). The 
coefficients are: 

The coupling constants k, A and y and the fields 0* and Xj dimensionless quantities. 
Np = 4 is the number of continuum four-component fermions. 

The scalar sector 5*^ of the action ( |2.1| ) has a global 0(2) symmetry. The action S is 
invariant under the vectorial U(N) transformations 

Xj ^ , xi ^ , e u{NpI2) (2.2) 

and the axial U(1)a transformations 

X - e'^^'-x , X Xe^"*^^ , <P ^ fi-^-A^ ^ ^* ^ e2..A^* , o^a e R . (2.3) 

The action (|2.1| ) contains two important limit cases, the model and the GN3 model. 
At ?/ = it is the 01 theory described by the purely scalar part Sb of ( |2.1| ). In the limit 



X ^ 00 the action Sb reduces to that of the XY3 spin model. At k = A = 0, the action S 
( p.l|) turns into the action of the chiral GN3 model in the auxiliary scalar field formulation. 
The full Yukawa model interpolates between both these models and the PM-FM critical 
sheet continuously connects the magnetic phase transition of the spin model with the chiral 
phase transition of the GN3 model. 



2.2 Symmetry breaking 



In order to get information about the breakdown of the continuous chiral symmetry in the Y3 
model we have computed the effective potential in 1-loop order for k > 0. For this purpose 
we start with the Euclidean continuum action with mo being the bare mass, go the bare 
scalar selfcoupling and yo the bare Yukawa coupling. The calculation is straightforward (see 
2^ ) and yields 
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where we have regularized the momentum integrals with a cut-off A. We have introduced 
the abbreviation a"^ = af + (t|, where the constants 0", (z = 1,2) can be identified with the 
expectation values (7^ = {(pi), (pi being the scalar fields in the continuum. These fields are 
related to the lattice scalar fields 0^ by 
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and the relations between the parameters are 
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All the values of (Xj which minimize Ves are possible candidates for the vacuum of the 
theory. We can find these minima by solving the equations dVes/dai = simultaneously 
for i = 1,2. One solution is cti = (T2 = 0. In the symmetric phase it is a minimum 
(5j5jVefr|(o,o) > 0), in the broken phase a maximum (9jc}jV^flf|(o,o) < 0) and a further solution 
exists. In this sense didiVes\(o,o) = is an implicit equation for the boundary between both 
phases of the theory. 

At fixed values of the parameters rriQ and Qq we can calculate the critical Yukawa coupling 



yc{nT'o,9o)- If we choose ml > we can always find a positive solution y^ of didiVcs\ 
which is 
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Equation ( p.7|) means that even for > 0, when the classical potential does not predict the 
symmetry breaking, a solution ydiriQ, go) exists. For all couplings yo with yo > yc the vacuum 
expectation value (ip) of the scalar field is nonzero and the chiral symmetry is broken. 

This computation of the 1-loop effective potential suggests that in the Y3 model, at 
sufficiently small n > 0, two phases of different symmetry exist, as indicated in fig. |l|. As 
usual, we call them paramagnetic (PM) for (0) = and ferromagnetic (FM) for (0) 7^ 0. 



2.3 The phase diagram 

Figure |I| displays a schematic phase diagram, including also some expectations for k < 0. The 
phases relevant for our purposes are PM and FM. In the PM phase both order parameters 
(0) and (xx) (fo^' 2/ > 0) ^-^e zero and fermions are massless. The lightest boson pair is 
degenerate. In the FM phase the vacuum expectation value of the scalar field, the chiral 
condensate and fermion mass are nonzero. 

To characterize the PM and FM phases numerically we have used the magnetization 
M = a/ + 4>1Yi V being the number of lattice points. A continuous phase 
transition is indicated by a singularity of the susceptibility 

X = F((M2)-(M)2). (2.8) 

For the numerical simulations we used the Hybrid Monte Carlo algorithm. The critical 
surface of the phase diagram has been found by localizing peaks of the susceptibihty on a 6'^- 
lattice. In table |1| the values of the coupling parameters at the maxima of the susceptibility 
are summarized. Of course, they give only an approximate position of the critical surface. 
In cases in which it was needed, the critical coupling in the thermodynamic limit has been 
determined by a finite size scaling analysis. 



Table 1: Peaks of the susceptibility in the Y3 model determined on a 6^ lattice 



y 




A 





1/6 








0.1730(6) 


0.01 





0.1818(5) 


0.03 





0.2008(10) 


0.10 





0.2265(15) 


0.30 





0.238(2) 


0.50 





0.249(2) 


1.00 





0.240(3) 


3.00 





0.218(2) 


00 


1.10(8) 





0.0 


1.25(8) 





0.5 


1.25(10) 





0.75 


1.28(8) 





1.0 


1.25(13) 





1.5 


1.25((13) 





2.0 


1.10(8) 





00 


0.3 


0.15(1) 





0.6 


0.12(1) 





0.80(8) 


0.08 






y 




A 


0.90(8) 


0.05 





1.42(7) 


-0.1 





0.0 


0.24(2) 


0.5 


0.3 


0.22(1) 


0.5 


0.5 


0.20(2) 


0.5 


0.95(13) 


0.1 


0.5 


1.10(10) 


0.6 


0.5 


1.30(13) 


0.0 


0.5 


0.0 


0.25(4) 


1.0 


0.3 


0.22(3) 


1.0 


1.05(15) 


0.06 


1.0 


1.30(25) 


0.0 


1.0 


0.0 


0.22(1) 


00 


0.3 


0.19(1) 


00 


0.6 


0.14(1) 


00 


0.80(8) 


0.1 


00 


1.0 


0.04(1) 


00 


1.10(8) 


0.001 


00 


1.47(8) 


-0.1 


00 



Our data strongly supports the expectation that for all positive values of y the condensate 
(xx) vanishes simultaneously with the magnetisation M. We have extracted the fermion 
mass from the fermionic momentum space propagator. The agreement with the tree level 
prediction amp = y{(t)) is quite good. In the FM phase we have also observed in the 0- 
propagator a massive particle, the cr-boson, and a massless particle, the Goldstone boson. 



These masses, as well as {xx)i ^I's not as convenient as the BCL cumulant for the study of 
the finite size behavior but can be used for a qualitative comparison of the physical content 
of the Y3 model in different parameter regions. 



2.4 Renormalizability properties 

Both the 0"^ theory and the full Yukawa model in 3D are perturbatively superrenormalizable. 
For the GN3 model this is different. The continuum 4-fermion coupling has negative mass 
dimension, and the corresponding interaction is therefore perturbatively nonrenormalizable. 
Nevertheless, it has been proved that the GN3 model is renormalizable in the l/A^iT'-expansion 

It has also been shown in the framework of l/A'^p-expansion that for weak scalar self- 
coupling A = 0{1/Np) the Gross-Neveu model and the full Yukawa model in 2 < < 4 



are equivalent field theories |T^, Near the nontrivial fixed point the kinetic term of 
the scalar field and the quartic scalar selfinteraction turn out to be irrelevant operators. 
However, in those works nothing beyond the range of validity of the l/A^iT'-expansion could 
be said. 



In ref. [O] the equivalence has been confirmed by analytic and numerical methods for 



the discrete chiral Z(2)-symmetry, still with A = 0{1/Nf)- We have extended that work to 
the U(l)-symmetric case and have investigated a wide range of parameters including infinite 
scalar selfcoupling. 



3 Finite size scaling theory 



3.1 The Binder method 



In order to examine the interplay of the universality classes associated with two different 
nontrivial fixed points in the Y3 model we have studied the finite size scaling behavior and 
tried to determine the critical exponents of the theory^ at several points of the critical surface. 
A very powerful method to do this is the Binder method of finite size scaling analysis of a 
cumulant [jl4|, |15l . 

It is sufficient to use scalar n-point functions even in the case of nonvanishing Yukawa cou- 
pling. We therefore follow refs. p!5| , and define the corresponding fourth-order cumulant 
Ul on a. cubic lattice of extent L: 
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where G^l^ and G^l' are 
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Xl,...,X4 l,J 



^The critical exponents ly, (3 and 7 are defined as follows {t is the reduced coupling): 

c~l^r^ M^t^ (fort>o), x~itr^- 



If both L and the correlation length ^ are sufficiently large then U l has the form 



Ur 



2 - 



/2 



(3.3) 



with analytic functions /2 and f^. Note that p.3| ) requires the validity of hyperscahng. 

At the critical value of the hopping parameter k, the correlation length diverges and 
all cumulants Ul\kc have the same value U* independent of the lattice size. This makes it 
possible to determine the infinite volume critical coupling as the common intersection point 
of Ul for different values of L. 

In the scaling limit Binder's cumulant has the form 



UL{t) = U 



L 



UitL-) 



\t\ < 1 



with t = 1 — Kc/ K. Let us consider a pair {bL,L) of lattice sizes with b > 1. From 
follows 

'dUtL 
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In order to obtain the derivative -gj^Uc o'ls calculates the function [7;,^ = g{U]}j numerically 
and near the critical point approximates (7 by a linear function determining its slope. 

Similar relations can easily be derived for the exponent 7 of the susceptibility x and the 
exponent /3 of the magnetisation M, 
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Using ( |3.6D one can calculate the ratios (ijv and 7/1/ from M and x determined on various 
lattice sizes (6L, V) exactly at k^. 

In the 0| theory the specific heat exponent a is negative. That means that the specific 
heat is a regular function of the reduced coupling and there is no relation similar to (|3.6| ) for 
it. 

To calculate the required quantities, we used a reweighting technique. By means of the 
original method due to Ferrenberg and Swendsen one can only interpolate operators 
which can be expressed as explicite functions of S. Therefore, like the authors of ref . |12[ , we 



used a variation of the method suggested in ref. . It can be regarded as the multihistogram 



method with bins of zero width. With that reweighting technique one can interpolate nearly 
arbitrary operators over a wide range of the coupling /?. For this purpose it is necessary to 
store the operator S which corresponds to the coupling /3 and the value of the operator for 
each configuration which has been generated during the simulation. 



3.2 Previous applications of the Binder method 

In the past the Binder method has been applied to a variety of interesting physical systems. 



In ref. the method has been generalized to 0{N) (f)'^ theories and in ref. the critical 
exponent v has been determined for the 0(4) invariant scalar 0"^ theory in 3D and 4D. A 
high precision measurement of u in the XY3 model has been done in ref. p8| . 



The method has also been apphed to models with interacting fermions [ll2|. Here a slight 



modification of the Binder method has been used to compute the critical exponents v and 
7/1/ in the 3-dimensional Gross- Neveu model with Z(2)-sjTiimetry. The found value of ^ 1 
is in good agreement with the prediction of the l/iVi?-expansion. 



In ref. [0 the critical behaviour of diluted Heisenberg ferromagnets with competing 
interactions has been investigated. The authors varied the concentration of spins and found 
two distinct universality classes which are separated by a crossover region. In this domain 
strong corrections to scaling appear, and Binder's method does not work well. Also evidence 
for a new, intermediate universality class has been found. 

3.3 Other methods to determine critical exponents 

For the 03 model we have also tried to compute critical indices by some other methods. 
Among these are the direct method, which makes use of the finite size scaling laws of physical 
quantities, and the scaling of the smallest Lee- Yang zero with the lattice size. 

On a finite lattice of extent L the susceptibility peaks at the value km{L) of the hopping 
parameter. If we increase the lattice size then km^L) approaches Kc according to 

km(L) -Kc oc L-^'". (3.7) 

Thus the measurement of km{L) for various lattice sizes L yields the critical exponent u 
by a corresponding fit. We have tried this method in the ^3 theory for different values of 
the scalar selfcoupling A. Our results were rather unsatisfactory because of their quite large 
statistical errors. For the same statistics we obtained more accurate values for u with the 
Binder method. 

Another possibility to determine u is to use the finite size scaling of the Lee- Yang-Fisher 
zeroes. By continuing the hopping parameter k, to complex values one finds that the partition 
function has zeroes in the complex plain. For finite lattices all the zeroes lie off the real axis. 
The zero kq with the smallest imaginary part scales like [^| 



Im(fi;o) ~ L 



1 



We have computed Im(Ko) on different lattices and extracted u from a double logarithmic 
plot. Our results are consistent with those obtained by the other methods but the statistical 
errors again turned out to be substantially larger than for the Binder method. 



4 Magnetic transition at vanishing Yukawa coupling 



4.1 The 03 model 



In the limit y = the action (2J.) describes free massless fermions and 0(2) invariant 03 
model with quartic selfcoupling. Besides our interest in the features of the 0| model as a 
limit of the Yukawa theory, here we have developed and tested the methods we wanted to 
apply to the more sophisticated and expensive fermionic model. The existence of a nontrivial 
fixed point and a finite nonvanishing value of the renormalized quartic selfcoupling Ar in 
the continuum limit make this model by itself very interesting from a field theoretic point 
of view, too. 

The phase diagram in the k-X plane, computed mainly on 6'^ lattices (see the y = 
entries in the table |I|), is displayed in fig. ^. The spectrum in the PM phase below the 
second order phase transition line contains two degenerate massive scalar particles. In the 



FM phase {k > Kc) the 0(2) symmetry is spontaneously broken and the hghtest particles in 
the spectrum are a massless Goldstone boson and a massive a boson. 




Figure 4: Phase diagramm of the </)| model. The A-axis has been rescaled: A* = 50A/(1 + 
5A). The circles show the maxima of the susceptibility on a 6^ lattice. The crosses show the 
positions of Kc determined with higher precision by means of finite size scaling methods. 



The renormalization group properties of the 03 model have been investigated e.g. in |30 
and are indicated in figs. ^ and ^ on the y = face of the phase diagram at k > 0. The model 
is superrenormalizable in weak coupling perturbation theory and its physics at infinitesimal 
scalar selfcoupling is dominated by the Gaussian fixed point (Gfp) at A = 0. At nonvanishing 
coupling A > the critical line Kc(A) is dominated by the IR-stable nontrivial WFfp. The 
investigations by means of e expansion or expansion of the 0{N) symmetric 03 suggest 
that the interaction term becomes irrelevant, and the only relevant term remains the kinetic 
one. This means that at A > only one parameter has to be tuned k — >■ Kc{X) in order to 
reach a continuum limit governed by the WFfp. Thus the same scaling behavior should be 
found when the critical line is approached at arbitrary A > 0. 



4.2 Results at A = oo and A = 0.5 

We have chosen A = oo and A = 0.5 and determined the renormalized coupling as well 
as some critical indices in runs in the k direction. A Monte Carlo determination of the 
renormalized quartic coupling \r has been done e.g. 



m 



3l[| for the Z{2) symmetric 



model. To our knowledge no analogous measurement exists for the 0(2) symmetric model. 



Following e.g. we define Xr in the symmetric phase as 

Xr = {LaniRfUL ■ 



(4.1) 



Here am^ is the mass of the a-boson extracted from the scalar propagator. To extrapolate 
to the continuum limit we varied the lattice size from L = 6 to 12 while keeping LaniR fixed 
to 4. 

At A = oo the renormalized scalar selfcoupling increases very slowly with the lattice size 
L. The linear extrapolation in 1/L to L = cx) suggests a value of A^ = 26 ± 4. At A = 0.5 



an extrapolation to L = cxd is less precise, suggesting Xr = 20 — 30. The agreement supports 
the expectation, that the model is dominated by the WFfp on the whole critical line A > 0. 
These results for Ar are also consistent with the expected theoretical value ||3^ . 

The most sensitive test for the appurtenance to the same universality class is the com- 
parison of critical exponents. Using the Binder method described in subsec. 3.1 we have 
determined the critical exponents (3/v and 7/1/. The method works very well at both A 
values in the whole range of lattice sizes used, 4^ — 16^. To illustrate this we show in fig. ^the 
determination of at A = 0.5. In fig. |^ the data for dUbL/dllL, used for the determination 
of u at the same A value, and the linear fit, are displayed. 
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Figure 5: The intersection point of the Binder cumulants in the (j)3, model on several lattices 
for A = 0.5 at Hc = 0.241(1). The lines were obtained by reweighting, the symbols are the 
measured points. 



A 




V 




7/z/ 


(") 


(5) 


(^) 


00 
0.5 


0.2275(10) 
0.241(1) 


0.673(19) 
0.687(19) 


0.51(3) 
0.56(5) 


2.03(6) 
1.91(6) 


-0.02(6) 
-0.06(6) 


5.2(4) 
4.5(3) 


0.02(6) 
0.12(10) 



Table 2: Critical exponents in the 4>\ model at A = 0.5 and A = 00. The exponents enclosed 
in brackets were calculated by using hyperscaling relations. 

We were able to determine i/ to a precision of about 3%, l3/u to about 9% and 7/z/ to 
3%. The results are summarized in table ^ They are consistent with the expectation that 
the two points A = 00 and A = 0.5 belong to the same universality class. The value oi (3/v 
is consistent with the one calculated with the hyperscaling relation j3/v = {d — 7/i/)/2. As 
hyperscaling seems to be fulfilled, we determined the exponents a, S and rj from the relations 

a = 2-ud, 6=P^, , = 2-d + 2^. 

d — 7/z/ z/ 



aUb,/aU,=f(b) A = 0.5 y=0.0 



linear fit: 

v= 0.687(19) 

c= 0.024(14) 



10° 



L =4 

N.,,„= 120000 



Figure 6: Linear fit to \n{dUf)L/dUL) used, according to eq. 3^, for the determination of 
the critical exponent v determined at A = 0.5. 



5 Chiral transition at vanishing scalar selfcoupling 



5.1 The GNg model 

At A = K = the scalar field (j) plays in the action (|2.1| ) the role of an auxiliary field. It can 
be integrated out thus obtaining a purely fermionic GN3 model with U(l) chiral symmetry, 



S — S p — 



y_ 
4 




beh.c. 



^x+bXx+bXx+b 



(5.1) 



In 3D this model is perturbatively non-renormalizable. However, it has been shown in 
that the GN3 model is renormalizable in the l/Np expansion. The /3-function has 
been calculated to 0{1/Nf) in |l^, |l9l and to 0{1/Np) in [^. The l/Np expansion reveals 
a nontrivial UV-stable fixed point where dynamical chiral symmetry breaking and fermion 
mass generation occur. The phase transition is of 2°^^ order and the order parameter is the 
chiral condensate (xx)- 

In ref. [|2^ one can find the critical exponent u to 0(1/ Np). In our case (A^^ = 4) 



1 + 



16 _ 8(3 76 + 277r2) 



+ 0(l/iV|) ~ 1 + 0.135 - 0.122 ~ 1.0(1) 



(5.2) 



The 0(1/ Np) term is identical with the results in |TP|, O]. The 0(1/^^1,) term is of the same 



order of magnitude, which suggests a rather large error on the value of v in ( |5.2| ). 

In the symmetric phase {y < fermions are massless. This region is dominated by the 
trivial Gaussian fixed point at ?/ = 0. 

By adding the kinetic scalar term to the bare GN3 action the scalar field (f) turns from 
an auxiliary field to a dynamical one. This restricted Yukawa model with A = 0, sometimes 
considered as a sufficient representation of the Fs model (e.g. in |jll|), is a natural extension 
of the parameter space of the GN3 model. We know that such a Yukawa model with vanishing 



scalar selfcoupling and Z{2) symmetry is renormalizable in l/Np expansion. As shown in 



10], this model has a nontrivial IR-stable fixed point where the kinetic term of the scalar 
field becomes irrelevant and the 4-fermion interaction term relevant. This fixed point is 
identical with the critical GN3 model. The IR-stable fixed point of this restricted Yukawa 
model corresponds to the UV-stable fixed point of the GN3 model 



lOl 11 . This can 



be understood from the renormalization group fiow in a larger parameter space, in the full 
Y3 model ( p.l|) (see figure 0). The fiow is suggested by the /^-functions obtained in the 
e-expansion ||^. The RG-fiow restricted to the GNa-line k = A = is consistent with the 
UV-stability of the nontrivial GN3 fixed point. 



5.2 Numerical results 

First we comment on the spectrum calculations. The fermion mass amp has been measured 
by fitting the momentum space fermion propagator, measured usually at four lattice mo- 
menta, to a free fermion ansatz. In the broken phase amp agrees very well with the tree 

2 

level relation amp = y{<P) = ^{xx}- 

For the measurement of the masses of the a-boson amo- and the Goldstone boson am^r we 
had to use an ansatz for the momentum space propagators from the one-loop renormalized 
perturbation theory |3^. In this case the previously fitted fermion mass is used to calculate 
the fermionic selfenergy which contributes to the renormalized bosonic propagators. This 
method delivers the renormalized Yukawa coupling and describes very well the form of the 
bosonic propagators which differ very much from the free ones. As expected, in the FM phase 
am.^ is very small and am^ increases with the distance from the critical point. In the PM 
phase both masses grow with the distance from the critical point and become degenerate. 

The scaling behavior has been investigated in two directions: In the GN3 case (k = 0) 
we varied y and determined the critical Yukawa coupling yc = 1.091(5) from the intersection 
point of the Binder cumulants on several lattices. From the finite size scaling behavior of 
the Binder cumulant at this value we determined the exponent u. Similarly the behavior 
of magnetization and susceptibility allowed us to determine P/u and 7/z/, respectively. The 
obtained results are collected in table |3|. 





Vc 


V 




7/z/ 


note 




0.000(2) 


1.091(5) 
1.09 


1.02(8) 
1.05(12) 


0.89(10) 
0.90(4) 


1.19(13) 
1.15(4) 


run in y (GN) 
run in k 



Table 3: The critical exponents z/, 7/z/ in the GN3 model and in the Yukawa model 
at vanishing A. 

By using the measured 7/z/ value one obtains from the hyperscaling relations j3/v = 
0.905(65). This is in good agreement with the measured value and supports the hyperscaling 
hypothesis. 

As a test of our methods and of the equivalence between the fixed points of the GN3 and 
the Yukawa model with vanishing A we measured the critical exponents in the latter model by 
approaching the critical point of the GN3 model along the k direction. Fig. demonstrates 
that the critical point obtained by the Binder method in this direction is identical with the 
GN3 one. The BCL cumulants intersect at t^c = 0.000(2) and yc = 1.09. As shown in fig. H, 
the values oi (3/v and 7/z/ are perfectly consistent with those obtained in the GN3 run (table 
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Figure 7: Determination of Kc in the Yukawa model at A = and y = 1.09. The intersection 
point of the BCL cumulants measured on different lattice sizes gives Kc = 0.000(2) which is 
in perfect agreement with the GN3 critical point. 
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1/ = 1.05(12) 



Figure 8: Determination of the critical exponent v in the Yukawa model at A = 0, y = 1.09 
and K = 0. 



We conclude that the Binder finite size scahng method is apphcable and gives consistent 
results in the Yukawa model at A = for a broad range of lattice sizes. The values of 
the critical exponents in the chiral GN3 model are the same as in the Yukawa model with 
A = 0. This confirms that the fixed points of these two models are the same. The exponents 
are consistent with the 1/Np predicted values (z/ 1) and significantly different from the 
exponents associated with the WFfp. This allows us to investigate the crossover effects 
between these universality classes numerically. 



6 Gross-Neveu-like behaviour for strong couplings 

The expansion predicts |10|, 11 1 that the Y3 model and the GN3 models are equivalent 
at least for weak scalar self coupling A. In order to test this hypothesis also for strong scalar 
couplings we have investigated the Y3 model with A = 00 at strong bare Yukawa coupling 
y = 1.1. This choice leads to — 0. 

The spectrum is similar to that of the GN3 model. We observe the generation of the 
fermion mass arriF which is related to a nonzero chiral condensate (xx)- Even for A = 00, 
where the l/A'^p-expansion is not applicable, we find that the prediction arriF ^ y{(j)) is 
fulfilled with good precision. Fig. ^ displays the dependence of the masses of both bosons on 
the hopping parameter k a.t y = 1.1 and X = 00. In accordance with the Goldstone theorem 
one massive a boson and one massless vr boson appear in the FM phase. The qualitative 
ft;-dependence of both masses in the vicinity of the critical point is the same as in the Y3 
model at A = 0. This is a first numerical hint for the physical equivalence of both cases. 
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Figure 9: The masses of the a- and vr-bosons as functions of k. As expected, in the broken 
phase the a boson is massive and the vr boson massless. 

In order to determine the universality class of the Y3 model at A = 00 and strong Yukawa 
coupling we have again determined the critical exponents (3/v and 7/1^. We have applied 
the Binder method aX y = 1.1 approaching the critical sheet in the k direction. The critical 
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Figure 10: An example (L = 6, 6 = 4) of the linear dependence of U^l on Ul at X = oo, 
y = 1.1 near the critical point. 



value Kc = 0.007(2) is given by the common intersection point of the cumulants Ul on 
different lattices sizes L. 

For this value of Kc we have computed the derivatives dUhL/dUL with L = Q and bL 
ranging from 8 to 24. Figure ^ shows as an example UbL as a function of Ul for 6 = 4. Near 
the critical point such functions are linear with good precision and the derivatives are thus 
easily determined. 

Using equation (|3.5|) we have obtained the critical exponent z/, 



0.89(6) 



(6.1) 



This value is a little bit smaller than the one obtained at A = 0, but both values are consistent 



within statistical errors. Figure |ri| shows the corresponding plot. We have also made various 
fits with different subsets of data points. The results are nearly unaffected if we leave out 
one or more data points in the fit. This shows that also for A = oo and y = 1.1 corrections 
to scaling are quite small and the Binder method works in a broad range of lattice sizes. 
We have further determined the ratios f3/i> and 7/z/, 



1 

V 



0.80(8), 
1.30(7) 



(6.2) 



Within statistical errors these exponents are consistent with our results in the GN3 model, 
too. They fulfill the corresponding hyperscaling relation with good precision. 

These numerical results lead us to the conclusion that the Gross- Neveu universality class 
extends over the whole range from A = to A = cxd provided the bare Yukawa coupling y is 
strong enough, y — 1. This confirms the conjecture that the GN3 model and the Y3 model 
are equivalent field theories even for A = cxd. 



aUbL/aUL=f(b) U(l)A = oo y=l.l 
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Figure 11: Log-log plot of the derivative dUbi/dUL at Kc at A = oo, y = 1.1. We show 
the linear fits to the first three,. . ., six points. The results of these fits are consistent values 
f3, . . . , of the exponent z^. 



7 Interplay of magnetic and chiral universality classes 

Both in the pure scalar 0| theory aX y = and in the GN3 model at A = the Binder 
method works in an exemplary way. Also in the Y3 model at A = 00, ?/ = 1.1 it provides 
satisfactory results. This is presumably due to the dominance of only one of the fixed points 
in these cases. They seem to be "pure" cases, without any interplay of universality classes. 
Now we describe what happens in the Y3 model when at A = cxo the Yukawa coupling y is 
decreased, and the XY3 model is approached. We made extensive simulations at y = 0.6 
and y = 0.3, approaching the critical sheet in the k direction. 



7.1 A = oo, y = 0.6 

For small lattice sizes, L = 6, 8, 10, the cumulants consistently cross in the interval k, = 
0.1460-0.1466 (fig. [l|a). Making the finite size analysis at k = 0.1463 we obtain u = 0.75(9), 
a value quite close to that of the XY3 model. 

However, when only large lattices L = 10, 12, 14, 16, 24 are considered, the crossing point 
is found in the interval k = 0.1430 — 0.1445. The situation is shown on a fine k scale in fig. 



TpD. For these lattices at k = 0.144 we find u = 0.99(23), a value consistent with the GN3 
model, but with a large error. 

We have made an analysis at = 0.144 including data on all lattices and choosing the 
basis L = 6. The u values have been determined for different groups of data, for the first 3, 
4, 5 and 6 points. As shown in fig. 0, when data on larger and larger lattices is included, u 
increases systematically from 0.71(11) for b = 8/6, 10/6 and 12/6 only, to 0.87(8) when all 
data is included. This is probably not a good way of analysis in such a complex situation 




Figure 12: Determination of Kc at A = oo and y = 0.6. On small lattices (a) the apparent 
critical k is 0.1463, but including also data on larger lattices and zooming into the rectan- 
gular region (b) suggests that the critical point is in the interval n = 0.1435 — 0.1445 of the 
intersection points of C//, for L > 10. 
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Figure 13: Similar to fig. 11, but now at A = oo, y = 0.6. The i/- values increase system- 
atically when data on larger and larger lattices is included in the linear fit. 



and the previous one made only on large lattices seems to be more reliable. We have made 
it in order to illustrate the systematic increase of the apparent v with lattice size. 

We interpret the above results as a hint that for sufficiently large lattices, L > 10, the 
xfp universality class finally shows up. It might be tempting to conjecture that the low value 
of obtained rather consistently for L < 10, is a signal for the nearby WFfp class. But, as 
the results at y = 0.3 indicate, this is questionable. 




Figure 14: An attempt to determine Kc at A = oo, y = 0.3 failed. Ul do not intersect in a 
single point even if only large lattices are considered. 



7.2 A = oo, y = 0.3 



As seen in fig. |1J, the cumulants obtained on lattices up to L = 32 show no tendency to cross 
at some unique point, even if smaller lattices are discarded. Also the dependence of Ubh on 
Ui^, shown in fig. |l^, is not linear, differing e.g. from A = oo,y = 1.1, seen in fig. |1^. A 
determination of u under these circumstances makes little sense, and one can only speculate 
that if lattices could be made still substantially larger, a simpler finite size behavior with 
the xfp exponents might be found. 

Remarkable is also the fact that the finite size behavior did not improve on small lat- 
tices. As in the y = 0.6 case, the cumulants on L < 10 lattices cross in a narrow interval 
He = 0.2045 — 0.2055.0 But including the L = 12 data spoils the consistency completely. 
Thus halving the distance from the XY3 model with respect to y = 0.6 did not increase the 
consistency of the finite size behaviour for smaller lattice sizes. This prevents us from inter- 
preting the low values of i' obtained on smaller lattices as a signal for the WFfp universality 
class. 

Attempts to incorporate some corrections to the leading finite size behaviour, as suggested 
in ref. are in our case not very helpful because simulations with dynamical fermions 

cannot yet produce data with the precision needed to deal with additional parameters. 
Thus, we conclude that at A = cxo and y = 0.3 the finite size behaviour is not under control. 



^The corresponding value of is 0.76 with errors difficult to estimate because of systematic uncer- 
tainties caused e.g. by a nonlinearity of the dependence of Uhh on C/l- 
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Figure 15: Example (L = 6, 6 = 4) of nonlinear dependence of C/fei on i/^ at A = oo, 
y = 0.3 near the critical point. 

Unfortunately, it would not have been easy to notice that without having data in a large 
range of lattice sizes. 

8 Summary and Conclusions 

We have studied the finite size behavior of the Y3 model with U(l) chiral symmetry along 
the 2 dimensional sheet of chiral phase transitions at various values of the Yukawa coupling 
y and of the scalar selfcoupling A. The aim was to investigate the influence and the interplay 
of the two nongaussian fixed points of the model for various values of the couplings. 

In the y = limit case, i.e. in the 0| model, the critical exponents associated with the 
Wilson-Fisher fixed point (WFfp) are clearly observed both at A = cxd and A = 0.5. The 
Binder method of finite size scahng analysis is applicable already on small lattices. Also the 
renormalized coupling values agree and are consistent with the theoretical prediction. The 
WFfp thus provides a rather complete description of the model at least for A > 0.5. 

For y > we find that the chiral fixed point (xfp) determines the finite size scaling in the 
vicinity of the chiral phase transition sheet as long as the Yukawa coupling is strong enough, 
y ~ 1. The independence on the value of the quartic coupling A confirms the expectation 
that the 1^3 model and its special case, the GN3 model, belong to the same universality class 
of the xfp. Also the fermion and boson masses at A = and A = 00 are very similar. For 
y — 1 the Binder method of finite size scaling analysis works consistently in a broad range 
of lattice sizes, in analogy to the pure (p\ theory. No substantial difference in the finite size 
behavior has been found between A = and A = 00. This implies that as long as y is large 
enough, the A0^ term does not influence the flnite size behavior of the Y3 model even on 
small lattices and the model is rather completely described by the xfp. 

When, at A = 00, y is decreased to y = 0.6, the finite size behavior cannot be analysed 
any more by the Binder method in the whole range of the lattice sizes we used. The behavior 
on small (L < 10) and large (L > 10) lattices is different. On the larger lattices the xfp 
seems still to dominate. On the smaller lattices the behavior looks consistent with the WFfp. 



But this does not necessarily mean that the WFfp already starts to show up: when a further 
step towards the 0| limit case is made, at A = oo and y — 0.3, the finite size behavior does 
not show increased resemblance to that fixed point. If applied in a narrow interval of lattice 
sizes, the Binder method might seem to be applicable but the results arc misleading. 

A numerical verification of the expectation that the Y3 model is equivalent to the GN3 
model is thus very difficult for y < 0.6. Our tentative conclusion is that the observed 
deviation from the simple finite size scaling signals an interplay of both universality classes, 
i.e. a crossover between them. This warns us that in the situation of intertwining phenomena 
the finite size behavior may be very complex. As we learned in the A = 00, y — 0.6 case, this 
fact is not easily detectable in a small range of lattice sizes, however. 
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